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ABSTRACT 


A  mathematical  model  of  a  branching  stochastic  process  utilizing 
generating  functions  is  presented,,  The  probability  distribution 
of  the  number  of  members  of  the  process  at  discrete  time  periods, 
zQ,  and  the  probability  of  extinction  is  discussed,,  When  there  is 
a  non-zero  probability  of  surviving  indefinitely,  the  normed  random 
variables,  zn/E(zn),  converge  with  probability  onej  the  cumulative 
distribution  of  this  random  variable  is  absolutely  continuous „  The 
time  until  extinction  and  the  total  number  of  members  of  the  process 
is  examined  when  the  probability  of  extinction  is  one.  The  distri¬ 
bution  of  the  zn  given  that  zn  is  not  zero  is  discussed  for  this  case. 
The  maximum  likelihood  estimates  for  the  probabilities  involved  in 
the  process  are  determined,,  An  example  is  given  of  a  branching  pro¬ 
cess  in  which  the  probabilities  are  dependent  on  time  and  a  solution 
is  found  using  Laplace  transform  methods. 
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NOMENCLATURE 


Symbol 

P(A) 
P(A  B) 
P(A,B) 


f(s) 

Pr 

w 

Pnr 


Gn(u) 

n(s) 

G(u) 

k(s) 

H(u) 

(s) 


Meaning 

the  probability  of  the  event  A 

the  conditional  probability  of  A  given  B 

the  joint  probability  of  A  and  B 

the  number  of  members  of  the  r™  generation 

the  number  of  offspring  of  the  i™  member  of  the 

n^h  generation 

the  expected  value  of 

the  variance  of  z. 

i 

generating  function  of  z-^ 

probability  of  r  members  in  the  first  generation 
generating  function  of  zn 

probability  of  r  members  in  the  n^h  generation 
probability  of  extinction 

z 

a  normed  random  variable  ( _ _ .) 

ri 

the  limit  of  wn  as  n  approaches  infinity 
the  cumulative  probability  function  of  wn 
the  characteristic  function  of  wn 
asymptotic  distribution  branching  from  f(s) 
probability  generating  function 
the  asymptotic  distribution  branching  from  k(s) 
the  characteristic  function  of  H(u) 
probability  that  zn  /  0 

probability  that  zn  equals  r  given  that  it  is  not 
equal  to  zero 
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yn  number  of  members  of  the  n^h  generation  given  that  it  is 

not  zero 

Fn(s)  probability  generating  function  of  yn 
N  number  of  generations  to  extinction 

(s)  moment  generating  function  of  N 

S  total  number  of  individuals  in  the  process 

g(s)  probability  generating  function  of  S 

zmr  number  of  individuals  in  the  generation  who  have 

exactly  r  descendents  in  the  (m  l)s^  generation 
Sn  total  number  of  individuals  in  the  first  n  generations 


vi 


INTRODUCTION 


The  present  paper  is  an  essay  on  the  theory  and  application  of 
discrete  branching  stochastic  processes.  Most  of  the  ssymptotic 
theory  was  extracted  from  the  work  of  Harris  [10.  who  derived  many 
of  the  relationships  for  this  class  of  stochastic  processes.  Let 
us  introduce  these  processes  by  means  of  a  historical  example. 

In  187b,  one  of  the  first  investigations  into  discrete  branching 
processes  was  made  by  the  Rev.  H.  W.  Watson,  Dr,  F,  Galton  had 
posed  the  problem  of  extinction  of  family  names  to  the  mathematical 
community  in  England  since  many  names,  prominent  in  history,  were 
no  longer  to  be  found.  The  problem  was  posed  as  follows:  Assume 
a  large  nation  with  N  adult  males,  each  with  a  different  surname. 

What  proportion  of  names  would  be  extinct  after  r  generations?  In 
the  rth  generation,  how  many  people  would  have  the  same  name?  Rev, 
Watson  [8]  derived  some  relationships  that  enabled  him  to  confute 
the  probabality  of  extinction  and  made  some  sample  calculations  to 
demonstrate  the  difficulty  of  computing  numerical  answers.  Notable 
among  these  relationships  was  the  functional  relation  of  the  proba¬ 
bility  generating  functions  involved  (section  III,  theorem  1), 

Branching  processes  received  some  attention  thereafter  and  the 
theory  was  applied  to  problems  of  gene  mutation  '7JS  population 
growth  [l2j  ,  cosmic  radiation  [l]  ,  ionization  fluctuations  [l8l , 
neutron  production  in  atomic  piles  [llj  ,  and  cascade  theory  [20]  , 

In  this  paper  we  will  review  the  theory  of  simple  discrete 
branching  processes  and  some  of  its  applications.  Let  us  consider 
the  original  problem  of  survival  of  family  names  and  state  the 
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basic  problem  of  branching  processes  in  that  framework.  Suppose 
an  immigrant  by  the  name  of  Door  comes  to  the  United  States  and 
there  is  no  one  else  here  having  that  surname.  What  is  the  prob¬ 
ability  that  the  name  survives  r  generations  and  how  many  men  will 
there  be  having  that  last  name? 

There  is  a  probability  that  Mr,  Door  will  either  not  marry  or 
not  have  any  male  children;  there  are  other  probabilities  associated 
with  the  number  of  sons  he  may  have.  We  will  denote  these  probabili¬ 
ties  by  pr,  i„e„,  the  probability  of  no  sons  will  be  pc,  of  r  sons, 
pr.  We  need  a  symbol  for  the  number  of  individuals  with  the  name 
of  Door  in  the  n^^'1  generation,  so  let  zn  be  the  number  of  individuals 
with  the  name  of  Door  in  the  n^h  generation.  Then  we  can  write 
P(z]_  ;  r)  =  pr  indicating  that  the  probability  of  r  individuals  with 
the  name  of  Door  in  the  first  generation  is  pr.  If  Mr,  Door  coes 
have  sons,  each  of  them  may  have  children  and  tt.e  process  is  continued. 
This  may  be  represented  as  a  tree  type  graph  with  an  initial  node 
representing  his  sons,  and  so  on,  A  method  of  analysis  of  the 
problem  based  on  graph  theory  was  proposed  by  Otter  [19]  ,  In  this 
paper  the  method  of  generating  functions  is  used. 

We  will  assume  that  the  number  of  sons  born  to  one  of  Mr, 

Door's  sons  follows  the  same  probability  law  as  that  of  Mr,  Door, 

In  other  words,  the  probability  that  Door,  jr„,  has  n  children  is 
pn.  Then,  if  is  the  number  of  offspring  of  the  i^h  individual 

in  the  n^h  generation,  we  have  P(zn+]_  =  r|an  =  m)  =  P(2r^  "  n) 0 

This  is  the  outline  of  the  problem,  the  notation,  and  the  assump¬ 
tions  made  in  this  paper.  Before  taking  up  the  mathematical  theory, 
some  other  problems  will  be  posed  in  this  framework. 
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In  atomic  piles,  the  number  of  free  neutrons  is  an  important 
factor.  If  we  start  with  one  neutron  and  a  fission  takes  place,  the 
probability  that  r  neutrons  are  released  is  a  nuclear  constant  which 
must  be  empirically  determined  for  the  amount  and  arrangement  of  active 
material.  The  probability  of  no  fission  (no  descendents)  is  the 
average  probability  of  leakage  or  absorption.  An  atomic  pile  in 
a  nuclear  reactor  is  subcritical  when  the  reactor  is  not  operating 
and  there  is  no  possibility  of  an  explosion,  i.e.,  the  probability 
that  the  number  of  fissions  increases  sufficiently  to  cause  an 
explosion  is  zero;  the  mean  of  the  probability  law  for  a  subcritical 
mass  is  less  than  one.  As  the  amount  of  fissionable  material  is 
increased,  the  probability  of  an  explosion  finally  becomes  greater 
than  zero.  When  this  occurs,  the  mass  is  termed  supercritical  and 
the  mean  of  the  probability  law  for  this  situation  is  found  to  be 
greater  than  one.  A  nuclear  reactor  is  termed  critical  Vhen  it  is  , 
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operating;  the  mean  of  its  probability  distribution  is  exactly  one. 

The  organization  of  the  paper  is  as  follows:  The  general  math¬ 
ematical  character  of  a  discrete  branching  process  is  presented  in 
section  I.  It  will  be  seen  that  the  mean  number  of  progeny  per 

individual  plays  an  important  roI,e  in  the  asymptotic  theory.  This 

> 

theory  is  presented  in  section  II  for  the  mean  of  the  probability 
distribution  greater  than  one  and  in  section  III  if  it  is  'not 
greater  than  one.  The  maximum  likelihood  estimation  of  the  distri« 
bution  of  progeny  per  individual  is  given  in  section  IV.  An  example 
of  an  ionization  cascade  is  presented  in  section  VQ  This  is  a  more 
general  branching  process  than  those  treated  in  the  model  used  in  the 
earlier  sections  and  its  treatment  suggest s  how  the  methods  of  Markov 
processes  may  be  applied. 
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GENERAL 


The  structure  of  branching  processes  was  presented  in  the 
introduction  where  several  assumptions  were  made  and  some  notation 
was  introduced*  It  is  appropriate  now  to  summarize  the  assumptions 
made  in  this  paper  and  to  extend  the  notation,. 

Consider  two  related  classes  of  random  variables: 

zn  the  number  of  M individual s,!  in  the  generation 
X^n)  the  number  of  "sons11  of  the  ith  individual  in  the 
nth  generation. 

They  are  related  by  the  equations: 

Zh  -  A  Xi  for  n=  1,2,3,.  o  . 

The  quantity  zQ  is  the  number  of  the  initial  group  who  start  the 
process;  zQ  is  assumed  to  be  one,  since  this  does  nob  affect  the 
generality  of  the  theory  as  is  pointed  out  at  the  end  of  this  section,, 
The  probability  distribution  of  the  first  generation  is  given  by  pr, 
i„e.,  the  probability  that  there  are  r  members  of  the  first  generation 
is  pr. 

The  following  assumptions  are  made  throughout  this  paper: 

(1)  Given  zn  is  equal  to  m,  the  random  variable,  zn*q  ,  is  the  sum 

of  m  independent  random  variables,  each  one  of  which  has  the  same 
distribution  as  zq  :  P *  =  r)  = 

(2)  The  variance  of  Zq  is  finite. 

(3)  1,  for  all  r  ;  otherwise  the  orocess  degenerates  to  a  simple 
arithmetic  calculation. 

(U)  pQ+  Pq  £  1,  otherwise  we  have  the  binomial  case  whose  treatment 
in  the  context  of  branching  processes  is  trivial. 
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The  following  notation  is  used  throughout: 


/4  =  Elz,)-S  rpr  , 


the  mean  of 


- 

e-'M?,)  = 

f!s)  =5  prSr 
r*o 

of  z^  (s  is  a  complex  variable) 


the  variance  of  Z]_ 

the  probability  generating  function 


Pnr  =  P<zn  =  r) 


OO 


of  zn. 


the  probability  generating  function 


From  the  assumptions  listed  above,  it  is  clear  that  : 
fQ(s)  =  s,  p^^  =  pr  ,  f 1 ( s)  and  f"(s)  are  continuous  on  the  set  of 
points  consisting  of  |s|<l  and  s=l. 

Several  authors  (see  [2],  [9],  jjll] )  have  examined  the  generating 
functions  and  have  found  some  basic  relationships.  Some  of  these  are 
presented  below. 

Theorem  1.  fn(s)  is  the  n^  functional  iterate  of  f(s)  : 


~  f  (Sj)  -  fn.|  (f  (s)J 


Proof :  From  the  initial  assumptions,  the  are  independent  for  all 

•  j  i~n  n,!1'1)  ,  \  ^n-i  M-i) 

i  and  n,  P  (X  t  =  K )  =  pK  Zn  =  2  X, 

u  1 


It  follows  that  P  (^n+1  is  the  coefficient  of  s^ 


in  the  power  series  expansion  of  (f(s))^  since  the  generating  function 
of  the  sum  of  j  mutually  independent  random  variables  is  the  product 
of  their  generating  functions.  See  Feller  [6,  p.  2!>0ff.] „ 
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Hence,  pA(s)  =  5  pj(f|Sl)  =  f(fU)| 

F,(s)  =  j>  R..[fls)lJ  =  ktflsl]  =FtfWsM  --  f  tfjlsfl 

*  j-0  J 


OO 


fn|s)  =  ^  pn.,Cf (sj]J  =  fM[f(s)]  =fn-.,|f  (.  .  (f/fisi))  . . .)) 

J-  0 

-  f  (Fn:,(s)J 

The  coefficient  of  sr  in  fn(s)  is  the  probability  that  there 

* 

are  r  individuals  in  the  generation,,  Thus  the  questions  posed 
by  Dr„  Galton  may  be  answered  in  the  framework  of  the  present  model,, 

The  constant  term,  pno,  is  the  probability  that  the  family  name  is 
extinct  by  the  n^h  generation. 

The  mean  and  the  variance  of  zu  can  be  determined  using  theorem  1„ 


Theorem  2, 


E  |Zn)=yUn  Varhn)  -  jul  *1 


f\C 


,  CO 

— :  Fn(«)  =  r Pnr  =  EUn) 


fnl ')  “^r(r-i)pnr  =  £(ZpJ~  £  (znj 


Therefore,  Year  ( 2  r )  —  ^  M  +  f n  (l)  —  , 

By  theorem  1  fn^j(s)  —  Fn  [f  (S)]  —  f  Cfn(S)J 


Taking  derivatives,  ~  [l"(s)jf  (S.)  —  f  L^n(s)]^n|S) 
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Fn",  is)  =  f"  C f  (s)3  [F‘(sO’l  +  tfcsijf  "(sj 
-  f"[f11lS)3Cf„‘(S)]i  +  f'Cfnls)]  C(S) 

Hence  Fnt|  t ) =  A  ^  dn)  ~ 


Solving  for  (1)  gives 

f"0)CfA<t»A-  f>f"W 

n  '  [  glllf  -  f‘l<> 

=  (6k-u+^Hu^-un) 

A*  -A 


sa/u.n(uri- 1) 

A^-a 


-  yU.”  tyU' 


Aa 


Thus  Var(zn)  =. 


g**"  (m"-i) 

Aa~A 


for  yx  ^  1 


For  the  case^t<  =  l  ,  L'Hopital's  Rule  can  be  used: 


d  n  “  I 


ei/U-l 


Am 


h 


-') 


6^  (2i o-n) 


A1' 


and  Var  (zn)  —  n  6"^* 


for 


si 


Higher  moments,  if  they  exist,  may  be  found  by  a  similar  pro¬ 
cess  utilizing  the  higher  derivatives  of  fnJ|.i(s)° 

* 

A  subject  of  great  interest  in  any  branching  process  is  the 
probability  of  extinction  of  the  species,-*  The  Fundamental  Theorem 
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of  branching  processes  tells  us  that  this  probability  is  a  fixed 

point  of  the  generating  function.  The  proof  given  follows  that  of 

Bharucha-Reid  [ 2]c  First,  we  note  that  if  "a"  is  the  probability  of 

extinction,  then  aslim  pnn=  lira  f-fo), 

n->oo  no  n-**5  n 


Theorem  3.  The  probability  of  extinction,  a,  is  the  smallest  pos¬ 
itive  number,  p,  such  that  : 

P:flph5  p,  pJ 

J-0  w 

Further,  a  a  1  if  and  only  if 

Proof :  We  see  ^t hat  a  :  0  if  and  only  if  pc  =  0  and  f(l)  =  1,  so 
that  0  —  Gt  -  1  as  it  should  be. 

To  show  that  a  =  f(a),  we  note  that  since  p.  =  fu(0), 


P-,,0  ■ 

so  a  =  lim  p  =  lim  fn(0)  =  lim  f(f  ,(0))  =  lim  f(p  ,  Q)  =  f(a), 
.  ri-»  °°  no  n-»<o  n-x  n-yoo  n-x,u 


Now  f(s)  is  a  power  series  in  s  with  positive  coefficients  and  hence 
strictly  increasing  for  real  arguments.  If  0<a<b,  then  f (a) ^ f (b) , 
Also,  since  prio^.O  for  all  n,  then  a2  0. 

Let  p20  be  such  that  p=  f(p),  It  follows  that 

p10  =  f|o)f  f(p)  =  p 

Assuming  pno^  p,  then  fn(0)£  p  and 

fmijO  ~  h)ti(o)  “.f  ‘  f(p)~  P 


So,  by  induction,  pno£p  for  all  n. 

Then  lim  pno£p  and  a- P« 
fwso 

It  remains  to  show  that  a  =  1  i'f  and  only  if  ^  ( 
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Assume  first  that  a  =  1  and  show  that  1 

Now,  it  may  be  noted  that : 

f'  (s)  Z  0,  f"(s)2  0,  so  f(s)  is  convex  in  (0,1) 


f(0)  =  po?0  ,  f(l)  =  1 

f(s)p»s  for  s  in  (0,1)  since  1  isthe  minimum  positive  number 
satisfying  s  =  f(s). 

Therefore,  1  -  f(s)<l~s  and. 

l-f(g) 


1-s 


is  bounded  by  1  and  is 
monotonically  increasing 
with  s  for  s  in  (0,1). 


Now  assume  Then,  since  f(s)  is  convex,  f'(s)  is  either 

constant  in  (0,1)  or  strictly  increasing  with  s.  By  assumption, 

^  1  so  PQ>0  and  f'(s)^  1  in  either  case  for  s  in  (0,1)„ 

By  the  mean  value  theorem,  P(s)  -  F  (Xl(l-S)  w here 


Then,  /-f(s)  =  f  '(v)(|-s) 

/-  f(s)  C  l-S 

F(S)>S  For  O^S^l 


and 


Therefore,  a  =  1. 

Using  theorem  3,  it  is  a  simple  matter  to  determine  the  limits 
of  the  pnr  for  r  /  0, 

Theorem  U.  A  discrete  branching  process  either  dies  out  or  becomes 

lim  p  _  =  0  for  r  i  0o 

n->cP  m 
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extremely  large,  i.e„. 


Proof :  Let  t  be  a  real  number  such  that  04  t<i  1  and  recall  the 


iteration  property  of  fn(t). 

Case  I.^U  41  Then  f(t)^t  (see  proof  of  theorem  3),  and 
f(t)Cf2(tK  .  .  .  ^fn(tK  .  .  .  <1. 

Therefore,  lim  f_(t)  =1=3, 

Case  II0/UVl  and  t^a  (pQ>0)„  Then  f(t)>t  since  f(s)  is 


the  limit  of  the  last  expression  gives? 

oO 

lim  f  (t)  r  lim  p  +  lim  ^  p _ tr  =  a  from  the  above  three  cases, 

n-*oo  n'  no 1  n->ao  ^ 

oo 

Therefore,  lim  >  p„,.tr  =  0,  Since  tne  series  converges  uniformly 
♦ 

in  an  interval  one  may  pass  to  the  limit  before  summing, 

oo 

Therefore  .•  lim  pnrtr  ®  0  and  lim  pnr  =  0  for  r  ^  0  by  the 

fT\  nr  nr 

uniqueness  of  power  series. 

Theorems  3  and  U  give  us  some  interesting  information.  If  we 
are  working  with  an  atomic  pile,  the  probability  is  one  that  it  will 
die  out  unless  it  is  sunercritical;  then  the  probability  that  it 
will  blow  up  is  1  -  a. 
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Lotka  [17]  found  that  the  probability  law  of  the  number  of  sons 
a  man  might  have  was  closely  approximated  by  a  type  of  geometric 
distribution,  namely,  pQ  =  0,h825,  p^  =  (0,212d)  (0,5893)k_1  k£Tl 

based  on  a  1920  census,  1  It  follows  that: 

V 

f  ( s)  =  Ooh825+  °-‘  21l6s~  U  ~  f'(s);  =  1,2622 

l-0.r893s  |  S=1 

Solving  f(s)  =  s,  it  is  found  that  a  s  0,6226,  Based  on  these  values 
the  probability  of  extinction  of  a  family  name  is  (0,6226) n  where 
n  is  the  number  of  male  members  of  the  family  at  the  present  time. 
This  formula  follows  from  the  fact  than.  if  there  are  r  members  of 
the  zeroth  generation  then  the  generating  function  for  the  first 
generation  is  (f(s))r.  Then  fn(s)  »  fn_]_( (f (s))r) 5  if  f(a)  =a,  then 
(f(a))r  s  ar  and  ar  is  the  probability  of  extinction. 
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MEAN  GREATER  THAN  ONE 


In  this  section  we  will  consider  the  case  where^  *  is  strictly 
greater  than  one  and  the  family  has  a  non-zero  probability  of 
surviving  indefinitely*,  Since  we  have  determined  that  zn  may 
become  extremely  large^  we  will  work  with  a  normed  random  variable* 


Let  wn  -  — 0  The  mean  and  variance  of  wn  are  easily  deter¬ 
mined* 


Elin)=l 


Treatment  of  the  convergence  of  wn  may  be  found  in  Harris  [lO] 
and  Bharucha-Reid  [2^)  0  We  will  prove  first,  the  convergence  in 
probability  and  then  the  convergence  with  probability  one. 

Theorem  If  the  random  variables,  wn,  converge  in  probability 

to  a  random  variable  w„ 

Proof,,  Let  n,  m,  be  integers  and  n>m<, 
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EU„|2m=r)=  rCF^lsr'Os 


(l)M 


/« 


-  TyU 

oo 

E  (2n  2,J  =  2  pm(.EPt>|*m»d 


-  2  P  r  V’"” 

r  wr  / 

I  yU"'"’  E<7m) 


5=' 


£  (Wn^m)  -  .  EiZftZw) 

r- 


_  — L  »-/oo  c  ly*.  \ 

~  JUL*"* /+  C  ?'«! 


E  (w£ ) 


c  [(u^-ww)1]  =  Elw’)-E(w^i 
..  <rA  /  j _ L) 

“  '  yiAm  jit*  I 

l"*>  Elu>»-w*f]  *  0 

m,  h-+<& 

Hence  the  random  variables  converge  in  mean  square  and  by  a  theorem 
of  Kolmogorov  ,  13,  p,3h,  Ij  ,  in  probability,, 

Convergence  with  probability  one  may  be  proven  with  the  aid  of 
Doob*s  f Uj  Martingale  convergence  theorem,, 
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A  discrete  parameter  martingale  is  a  stochastic  process  such  that; 


a)  E  (}yj  )  £  co  for  all  n 

b>  E  ^mll  7n’  yn"1,  *  *  °  >  yl^  s  yn 

The  theorem  may  be  stated  as  follows.  Let  (yn,  n20)  be  a  discrete 
parameter  martingale.  Then  E( |yQl  ) £  E(|yfl  ) £  E(  )~  »  *  » 

If  lim  E(  yn(  )  =  k  £.  °0  f  then  lm  yn  s  y*  exixts  with  probability  one 
and  E(|y*|  )  £  k. 

Considering  our  random  variable  wn  ,  E(wn)  =  1  for  all  n, 


'  yUn 


-  W„ 


Hence  the  (wn)  form  a  discrete  parameter  martingale  and  we  have 
Theorem  Ha,  If  ju> 1,  the  random  variables,  wn  ,  converge  with  probability 
one  to  a  random  variable  w. 

Since  wn  converges  to  w,  we  may  discuss  the  probability  distri¬ 
bution  of  w„ 


«3 

Let  Gn(u)  =  P(wn*u),0n(s)  =  E(  ewns)  «  Joesu  dGn(u), 

G(u)  2  P(wi£  u) ,  $(s)  2  j  esu  dG(u),  The  distribution  G(u)  can  be 
called  the  asymptotic  distribution  branching  from  f(s).  There  is  a 
very  interesting  functional  relation  between  the  moment  generating 
function  (JP(s)  of  w  and  the  probability  generating  function  f(s)  as 
given  by  the  following  theorem. 
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Theorem  $„  (ft  (yU  S)  =  fUftiS)’ 


Proof : 


ft  in  ( s )  =  E  (e^J  -  t  :=- 


X. 

J-0 


=  fn(e^j 

-  rCCRn(s)3 


6/}(s)  r  6  (s)  , there fore, 

M-»co 

A's.  ^ is>  ’jiz,  cv^  m 

=  J0VudSM  =  a>ls) 

Upon  taking  the  limit  of  iPnil  (yUS)~  '  [  (SO  the  desired  result 
is  obtained,, 

It  also  follows  that  the  k1-*1  moment  of  w  exists  and  is  finite 
if  and  only  if  the  k^h  moment  of  exists  and  is  finite,, 

Since  Gn(0)  s  P(wni  0)  =  pno  ,  then  lim  Gn(0)  "  a  and  P(w  =  0)sa 
The  following  theorem  will  be  used  later „ 

A. 

Theorem  6„  Let  G-^(u),  G2(u)  be  distributions  having  equal  first 
moments  and  finite  second  moments  such  that  their  characteristic 
functions  (fi  ( it)  and  (^2  (it)  satisfy  p  f  (fp  (^)J  r  =  1,2 

Then  G-^(u)  s  GgCu) 
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Proof : 


K-l 


if  tit)  =  |  +  "d  (it  A  o(iK) 

V-i  V!  ' 


for  ft  small  by  a  special  form  of  McLaurin's  theorem  for  small 
values  of  !t|  .  See [Cramer  3,  p»27]  .  Then, 

(it)  '  f  Q (t^J  „ ^ 


l U>aUt) I  =  (  1  +  "v-  «l 


r.  i 

a.  (3)  nsl 

2  ^f  tu)'1  t  O  (t3)  where  If,  '"lit)  - 


Itf.kt)-  cfe(it)|«  lhuU,uV 

~  +  o(t)l 

Therefore,  J  iff  lit)  -  (l  i  )  j  ~  £  (tj  Ujl^eirg 

Then  J  tP,  UiyU.)  -  *  |  f  Ctf|  Mff-  by  theorem  5 

£  yU.  |  Lpj  lit) -(A  1^)1  since when  |s|£j 

Hence  j  i‘lyU'L^  (yxi)  (  £  y(/t|£V^££)| 

and  yu|^(yl/tt)J^  |(3(0|  implying  that  yOtl^(i) I  -  I  |#(oJ| 

But yx.  >1  so  (t)  must  be  identically  zero„ 

Therefore  j  $,  ( U)  “  4>AU£)|  -O 

and  ( it)  s,  lit) 

Finally  Gq_(u)  ”  G2(u)  by  the  uniqueness  of  moment  generating  functions,, 

In  order  to  determine  some  properties  of  G(u) ,  it  will  be  con¬ 
venient  to  work  with  a  closely  related  function  H(u)  which  is  defined 

as  the  asymptotic  distribution  branching  from  k(s)  where 
j^(S) -  -fi 


!  -  Ca 


=16. 


Mp 


To  demonstrate  that  k(s)  is  a  generating  function  use  the 


fact  that  yU>l  implies  f(a)  ~  1.  Then 

'f'l  o)s  =  o 

!<(/]  -  14^-—  =  i 


Clearly  k(s)  is  monotonically  increasiiig  since^,  if  t^s,  then 

f  [111— LtgQ  vr  r[s(i-Q)-tOl  -q 

/-a  /~ci 


Therefore  k(s)  is  a  probability  generating  function,, 

It  is  interesting  to  note  the  following  properties* 


i<)  ~  hi!  '  -aj  _  f 1  fs  (t-al  f  a] 

'  /-a 


K’li)  -  f'(i)  yw. 

!s)=  f"Csti-a|+a](l-a) 

K"  1 1)  -  (sA+/tt.3-/u.)(i-a) 

K"  (1)  +K’0)- ( !)la  =  S^U-a)- a  l/Jyti) 


K(s!  = 


1 

j-o 


00 


OO. 


^0 


The  coefficient  of  s^  is  (l-a)^~^(  D  x  P  K&  •*'  ,  .  ,  +  is  )  P  G.'  + 

l/K  >k*t  '  A  f  lx+n 

and  C|j  =  f)  +  olp^a.  +3  f»aA  t  .  .  .  =  f’tal 

Let  l^(S)  ^  f  Sj£*sd  L'  uJ  o  By  a  line  of  argument  similar  to 

the  proof  of  theorem  5,  it  can  be  shown  that  ^  (ytt  S)  =  K  C  lp(S)]. 
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To  demonstrate  that  k(s)  is  a  generating  function  use  the 


fact  that  JU>1  implies  f(a)  z  a^  1,  Then 

k  to)  -  .  o 

t  — '  C4. 

Kh)  -  ■  =  i 

I  -  (X 


Kit ) 


Clearly  k(s)  is  monotonically  increasing  since, 

f  n  (i-o)fa]-a 


/-a 


[[slbaHd]  —  a 
/-a. 


if  tk  s,  then 

“  k'u) 


Therefore  k(s)  is  a  probability  generating  function 

It  is  interesting  to  note  the  following  properties,, 


ft1  (s)  -  -n •'  - ’ .  -  f  sfs(f~aH<2] 

i  ^  Ol 

!<'(/)-  "  '  (  ')  -yM. 

K''  =  f“[  slt-aH-a]  ll-al 

K"  til  •  (6* 

K"  ( i )  +  K’OI  -  [K'for  -  S2II-a  )-a^yut 


K(s) 


oo 


r=o 


r 


0 >0 


The  coefficient  of  sk  is  (l-a)k“1(  p  P  K&  +  ,  .  ,  t  j  l/  ]  P,  .  »  . * 

.«K  1  A  >  iKtft 

and  qt  =  p(  t  ap^o.  *  3  *  .  .  -  =  f  '{ft). 

Let  y/(S)  =  f  U  tu)  c  By  a  line  of  argument  similar  to 

the  proof  of  theorem  5,  it  can  be  shown  that  {jU~5>)  =  K  Cl S)J , 
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We  now  proceed  to  determine  a  representation  for  the  random 


variable  we  First  we  determine  the  relationship  between  W  (s)  and 

(P(s). 

Theorem  7. 


^(sj=  je£Srr£~<l— 


Proof:  Let  %($)&  (Mz&lSJ  ~  O. 

/  -CL 

Then  fj/,'  (S)  a  iOlLt^lLl^L .  -  W!f  <]-&.)  s] 

!~<x  *  * 


and  —  cP  ( O)  -yK.  o  Therefore  and  have  equal 


first  moments.  Since 


s  ip  [0-a)s]  (t-&)|  <  ^ 


and  Jr  have  finite  second  morents. 


iu,iits)=  ipcd-aisMi-Q-  -  f[ 

ri  7  y-a  /-a 


a 


-  rC(i-fl) W (s)-f-g]  ~  a 

/-a 

tp>3  (jus)  ^  \<i%  Is)] 

Thei’efore  y  and  $  satisfy  the  same  functional  relationship  with  k(s), 
Therefore  by  theorem  6,  ^{l}  =  Ip  ($}, 

Now  we  can  determine  the  relationship  between  H(u)  and  G(u). 

/  y  \ 

Theorem  8.  ft  { u)  -  6  [1-a  }  ~  0- 


/-  a 

Proof  g  From  theorem  7,  ^  (S)  “  (,PO?°Qj53°Q- 


-a 


By  definition  of  generating  functions  {jjj  ( S*  s  j  d  H(u) 

and  -_L_r  «*  Jg  (w)  _a  1  =  JfMl&JSlza 

/ “ o.  L  Jo  J  ’-ci 


Hence 


r»o 


<?SUJWM  =  yfa  [  £y“'a)s  JGH  -  a] 
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S~U-a)e?uJHM  =  J“e^"a)sJG(uj-a 

sj  0 


Let  f(t»J  be  defined  as  follows:  £  j  U )  -•  0  r  Cf  U  £  O 

l(fj)  “  Cl  C0y  70 

Then  O.  “  )  d  £(u)  “  f  _  J5  dc(d  ‘ 

J  Qm  J 

Therefore  a  is  the  Laplace  transform  of  £  ('Jl. 

Hence  j  £S°  (h<X)  d  H M  “  ,  S  ^ d  G  ("fT^}  ~  £  $  Cf'1’ 

and  CesuJit>-a)HM]  ~  Ces"J[S(-,^s)‘ €(!■•)] 

'JO*  '& 

Therefore  (l-OiHN  “  6  (  j  “  £M 

and  W(t»)  a  j  6  (7^)- a]  tor  vrO 

Clearly  H(u)  a  0  and  H(«* )  =  1  as  it  should<> 

To  put  the  relationship  of  H(u)  and  G(u)  in  terms  of  random 
variables,  we  prove  the  following  theorem0 

Theorem  9 »  The  random  variable  w  is  distributed  as  the  product  of 
two  indepei  dent  random  variables,  wQ  and  w',  where  P(w0  ®  0)  ~  a, 

■j 

P(w,_  s  - - )  s  l-a  and  w#  has  t»he  distribution  H(u)„ 

°  l~a 


Proof : 


E(eHUJ‘sj  =  a  +  <i*at£%  * 


-  ct  4  U-a)^  e H  (vU-d)\ 

© 

*  f  e'Jc/(SW 

"'O 

s  CP  ( S) 
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is  given  with  some 


So,  P(  ui(i-Cl)  f  U  |i^>o)  -  H(u). 

For  later  use  a  theorem  due  to  Harris  [  lo] 
remarks.  Let  J  =  lo^  L.)  =  log^^-L^).  If  ^  =  0,  let  7-  oo  . 

Theorem  10.  If  ^cco,  Re(s)^0  ,  s  /  0,  then  ^(s)  =  ■  -y  *f  M0(s) 

9  ^  f  s 

s)  =  M(s);  M  (s)  9  0  Mg-J  .  IS l-»°° 

a  S I  / 

M(-s)  =  M(  s) . 

Remarks:  (a)  Under  the  conditions  of  the  theorem  M(s)  is  real  and 
poaitive  when  s  is  real  and  negative. 

(b)  If  p  ( j  <:  c>o ,  the  r^’  derivative  of  li^(s)  satisfies 


(c)  If  If -z  oo  ,  U^(s)  and  as  many  derivatives  as  exist  approach 
zero  exponentially  as  |s|-»do. 

We  will  now  determine  some  of  the  properties  of  G(u) ;  however, 
it  will  be  convenient  to  orove  theorems  for  H(u)  and  then  interpret 
the  results  in  terms  of  G(u) .  Let  h(u)  =  h'(u)0 
Theorem  11..  H(u)  is  absolutely  continuous . 

Proof :  Cramer  [3]  breaks  up  an  arbitrary  cumulative  distribution 
function  F(x)  into  three  pieces: 

F{X)  ~  Q|  F,(x)  i  aA  FaM  4  F3  lx) 

where  &  f  -j*  a  3  =  1  0.^  Z  0  j  l  -  \j  2,  3 

F^(x)  is  absolutely  continuous. 

F2(x)  is  a  step-function  and  is  equal  to  the  sum  of  the  saltuses 
of  F(x)  at  all  discontinuities  which  are  less  than  or 
equal  to  x„ 

F^(x),  the  "singular"  component,  is  a  continuous  function  having, 

almost  everywhere,  a  derivative  equal  to  zero. 
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where  M(s)  is  continuous  for  s  f  0,  M (JJ. 


QO  ,  ^  y 

st)  ~  p_  Q  w  is  the  sum  of  an  absosutely  convergent 


h  =  0 


to 


trigonometric  series  and  is  thus  an  almost  periodic  function  which 
comes  as  close  to  ag  as  we  please  for  arbitrarily  large  values  of  t„ 

2 

m  -itu 


By  theorem  10,  lim  IIJ  (t)  =  0o  Therefore,  for  H(u),  a9  “  0  and  H(u) 
is  continuous a 

Let  hm(u)  =  “  f"  eLtv  UXit) 


Let 


l  m 

-ciu 


dv  =  e'^Jt 

I  'itu 

v  =  -  A  e 


u=  (j/Ut) 
du=  J[  V(U)] 


L  (u)  = 


n 


m 


ollT  L  U 


eitu  iput) 


DO 


+ 


_  J 

c2?n  u 


■  itu 


vr 

j  ,M  i  v  ni,  > 


<j  (put) 

(J  t 


Jt 


(A)  =  “  377 IU  | 


,  mo  •  ( t  U 

__e _ 

Jt 


dl|M 


Examine  each  term  of  (A)  as  m  oo 


First  '\]/ltnt)e 


-  imu 


M  Urn) 

l  ml * 


t  M 


i  ll 

0(s)je 


m 


i  Oi'jJ. 

I  ml  j 


by  theorem  10 . 


As  we  pass  to  the  limit,  this  term  is  a  function  of  continuous 
functions]  hence,  it  must  be  continuous.,  Similarly  the  second  term 


leads  to  the  expression  ip 


-i  mu 


I  ml  M(ur>)  t  0 


h»} 


i 


which  converges  to  a  continuous  function0 

.  rr  ^itu 

The 

» *>i 


e  last  term  in  (A)  is  i  y - —  jj  J* 

2-mu  d  i 

Since  |  J  p/y) dx  [  £  J  JpW  6x  ,  it  follows  that 

—  -itu  I  nil  .1  , 


r M 


(j  (put) 


J  {put} 

£rLU  at 


J  t 


*  K  m 0\ 


6  ><  0 


( i  Jf ) 


by  corollary  b  to  theorem  10 
for  some  constant  k. 


Hence  this  term  approaches  zero  as  m-^^o  0  Therefore,  the  continuous 
functions  hm(u)  converge  uniformly  in  some  interval,  u^  ^  u  ^  Ug, 
to  a  continuous  function  h(u)„ 


)  "  H(l/J  =  lini  r 

no-?co 


►n  ~tt  u, 

s  -  e  -  w)  d± 


Q.  TT  it 


'  .'/TO  —  * 

ATT 


^  m  r 

J-m  Jl 


-t£ 


u 


di/dt 


/, 


‘rri 


-Ms  i . .  r  r»  -  l tv 


r a  r 

Ju,  *tt  J„ 


LplitJ  Jt 


since  b^(u)  converges  uniformly 


JUf  3-TT  J.jy, 

”  ri'^U  jUi  h^(u)du 

to  h(u)  the  limit  may  be  taken  inside  the  integral  and 

HM-HM  =  ]■“•  AMJu 

Since  H(u)  is  continuous  at  zero  and  we  have  shown  'bsolute  con¬ 
tinuity  elsewhere,  it  is  absolutely  continuous  on  an  interval  in¬ 
cluding  the  point  zero. 

In  case  £  ,Z,  and  r2  JT+K-I ,  integration  by  parts  of  (A)  and 

corollary  b  to  theorem  10  shows  that  the  first  r  derivatives  of 

h(u)  are  continuous  if  u  /  0o  -  The  integral  expression  for  h(u) 

(T\ 

in  terras  of  ^(it)  shows  that  (fpr+i  implies  hv  y(u)  is  continuous 
at  zero.  This  implies  that  G(u)  ”  a  -f  f0  for  u?0;  g(v) 

is  continuous  for  v  £  0.  If  then  g(r)(u)  is  continuous  for  u^O 

provided  r<jr+/<-/  end  is  continuous  for  u  ®  0  provided  r<;  y- 1  „ 
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Corollary  to  theorem  11,  s  P(zn  -  r)->0  uniformly  in  r?  r2  1 

as  n^yoo, 


Since  Gn(u)->G(u)  uniformly  for  u>0,  G(u)  is  uniformly  continuous 
for  u>0  and  G(u)  is  right  continuous  at  u  s  0>  the  corollary  is 
provenQ 


It  is  interesting  to  note  that  Yaglom  [22]  treated  this  subject 
in  a  manner  similar  to  Harris .  Yaglom  defined 

S  Pro 

Vr  r  =°nw« 

M  “  P(yn  *ylyn>o) 
y(y)  =  hm  Hnlyi 

^kt)=f  eltyJ  HM  and  proved  the  following  theorem0 

—  OO 

Theorem  1$„  If  jU>i  ,  the  variance  of  z ^  is  not  zero,  and  f"(l)^co 

then  fttyUt}  0-flJ  +  &J  5  (l-ClJ  id 

■  '  In  this  paper  we  followed  the  general  line  of  Harris  argument 

and  defined  k(s)  such  that  ty(lis)  s  ~ — ~ — — —  , 

ff  dt 

This  led  us  bo  the  conclusion  that  II(u)  was  the  distribution  function 

2» 

for  (l-a)w  which  is  the  limit  of  (l-pno)-$-  0 
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Theorem  17  „  IfyW.£|  and  f  <  c*3  then 


!>  m 

p* 

-  p* 

n-^oa 

lr.tr 

r 

1  >  m 

fvts) 

~  F(s 

r--l 


and  fCf  'S))-yii  P/s)  +  {fyu/  jsl^) 


Fill- 1  r‘(l)*K'J 


Proof : 


p(z„=r  7n*o is" 


oo 


Pll^lLZaig)^ 
P(  z„«)  S 


by  definition  of  conditional 
probability 


-  ”  .r 

fr;  P(Z„  ■#£)) 


since  the  summation  is  over  r  0 


-  V  E&irj  s  r 

_  _£,  fs)  -  Uc) 

~  -  fnlo)~ 


(1) 


L/jItJ 


The  following  theorem  may  be  found  in  Konig  ^lU'o  If  5>(s)  is  analytic 


in  a  domain,  such  that  0(sQ)  =  sQ  and  q\sq)  “  a^l  for  some 
sQ  in  R,  then,  in  this  domain  —  “ 


a 


converges  uniformly  to 


a  continuous 


function  D(s)  such  that  D  -  Q  <J 


f(s)  is  analytic  in  ,  s  <  1.  Since  f(l)  a  1  and  f'(l)<  1,  (/k<Cl), 
the  sQ  of  Konig's  theorem  is  one0  We  can  conclude  that  f ^ (S '  -  1 

An 

converges  uniformly  to  a  continuous  function  D(s)  such  that 


(2)  D(f(s))  =/£D(s),  D(l)  =  0  ,  D'(l)  s  1„ 


(3)  F(s)  =  1  +  D(s)  K"^  by  Konig's  theorem,  since  convergence  of 

■£L-:j^ — ; - is  uniform,  the  limit  of  exists0 

J ^4.  yU- 

From  (3),  then,  DO)  —  (  F  S)-  )  K 
From  (2),  ~  DO) 

=  julKZ  F(s)-I3 

FCffc}3  =  I  i  by  (3) 

K 

-  ■  +  /  }< 

“  !  +  yU  FI  S)  ~/A- 


And,  finally,  “Cc(S>j  =  F  (s)  f  ( I'yuj 
From  (3),  F(j )  -  I  t  DO)  K~’ 

—  |  since  Q  (  | )  =  o  by (2) 
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One  of  the  more  interesting  questions  in  the  subject  of  branching 


processes  is  the  determination  of  the  number  of  generations  to 
extinction  when  yU<l» 

Let  M  -  mm  Ln  ;  zh*(  ~0j 
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Theorem  18,  The  probabilities  bn  satisfy  the  recursion  formula,. 
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Theoren  196  The  moment  generating  function  of  N  is  given  by 


&(%\  ~  (  +  ~  1  j  S|  (  s) 
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So  @(s)  is  determined  once  £^(s)  is  known.  If  we  let  z  -  e5, 
£)^(s)  ^"S  a  power  series  whose  coefficients  are  successive  iterates 
of  f*(b)  =  1  -  f ( 1-b)  since  bn+1  =  f*(bn)  =  fJAl(b0),  where  f*(0)  s  0 
f*  (o)  io  (see  Fatou  L  8  J  . )  A  function  of  this  sort  is  meromor- 

phic  with  poles  at  s  s  -  n  logju  ,  n  *  1,2,.  .  .  So 

®,(s)  iyute*  '^>es 

converges  everywhere  except  at  the  poles,  (see  Lattes  [29]  ) 

OQ 

X(s)  -  x^s^  is  determined  by  X(  a  s)  "  f*(X(s) )  ,x'  (1)  -  Xy(l)s  1 
y0  is  found  by  X(yD)  s  b0  s  1  »  pc  or  by  using  inverse  functions 
X=>l(f^'(s))  s  y6dX“f(s)  and  determining  X"^(b0)o 

If  our  process  is  of  the  type  in  which  the  "parent"  does  not 
die,  then  the  total  number  of  individuals  is  of  interest.  If  ja  £  1, 
the  probability  that  this  number  is  finite  is  equal  to  one  by 
theorem  ha. 

Let  ?«/,  then  p/S<cP)^  '  for  I 

■*=°  no 

Let  O  -P(s~r)j  for  *S  <T. 


rsc 


Theorem  20.  Ifyui  1,  then  g(s)  =  s  f(g(s)). 

Proof?  Let  2*  for  some  n  with  probability  one 

when  //  £  1  by  theorem  Ua.0  Therefore,  Sn  converges  to  some  random 
variable  S  with  probability  oneG  The  random  variable,  represents 
the  total  number  of  individuals  in  the  generations  from  the  zero^*1  to 


the  generation*,  Let  $ where  Vy 


.  r  ? 


since 


Let  £„{S' 
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where  =  c  and  -  P ( r-  r)  . 


The  event  s  n  can  be  divided  into  n  mutually  exclusive  events 
with  probabilities : 

P(z1  “  i)P(n-i  individuals  in  the  following  k-1  generations) 


for  is  1,2, 
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So  @(s)  is  determined  once  £^(s)  is  known,  If  we  let  z  ~  es, 
g-^(s)  is  a  power  series  whose  coefficients  are  successive  iterates 
of  f*(b)  =  1  -  f(l-b)  since  bn4l  S  f*(bn)  =  f*Vil(b0)s  where  f*(0)  =  0 
**  (0)  ~p*<  10  (see  Fatou  £8j»)  A  function  of  this  sort  is  meromor- 
phic  with  poles  at  s  =  -  n  logyU.  ,  ns  l#2r  0  „  So 

S/s,=  at  +  .. 

i -/<«*  ye‘  '->7e 

converges  everywhere  except  at  the  poles0  (see  Lattes  ■  29]  ) 

oO 

X(s)  x^si  is  determined  by  X(^s)  ~  f*(X(s)),x'  (l)s  X'(l)s  lc 

i*T 

yQ  is  found  by  X(yQ)  r  bQ  s  1  -  pQ  or  by  using  inverse  functions 
X~l(f#(s))  s  yWXr^-(s)  and  determining  X"^(bc)„ 

If  our  process  is  of  the  type  in  which  the  "parent"  does  not 
die,  then  the  total  number  of  individuals  is  of  interest.  If  u  £  1, 
the  probability  that  this  number  is  finite  is  equal  to  one  by 
theorem  Ua„ 

Let  $s]>~  2V>;  then  p{S<(cO)  ~  1  for  ^ j 

CO  ^ 

Let  qr  =  P(S~r.'^  for  15  * 

f  ~c 

Theorem  20 „  If  ft  £  1,  then  g(s)  £  s  f(g(s)), 

tn 

Proof:  Let  '  >  1*  2»  -  0  for  some  n  with  probability  one 


L~o 


when  //  £  1  by  theorem  Ua,  Therefore,  Sn  converges  to  some  random 
variable  S  with  probability  one.  The  random  variable,  Sn,  represents 

the  total  number  of  individuals  in  the  generations  from  the  zero^n  to 

* 

the  generation*,  Let  where  ?  since  7b~ 

Let  $w  ($  =  /  ,  where  =  o  and  2y>  =  Pi Vr-r). 

n=o 

The  event  y^  s  n  can  be  divided  into  n  mutually  exclusive  events 
with  probabilities: 

P(zx  s  i)P(n-i  individuals  in  the  following  k»l  generations) 


for  ±s  1,2,0 
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Starting  with  one  individual  at  time  zero,  the  probability  of  pro¬ 
ducing  n  individuals  in  the  next  k  generations  is  the  coefficient 
of  sn  in  G^(s)0 

IT- 

Therefore,  q.  a  .  ?(y->  *  i)P(n-i  in  the  next  k-1  generations) 

Kj  T1  J- 

x  . 

=  '>  p^( coefficient  of  s1*”1  in  (^^(s))1) 

"  \  p^(coefficient  of  sn  in  (sG^.^^Cs))1’-) 
s  coefficient  of  sn  in  (  Pi(s^k-l(s) )^) 

"  coefficient  of  sn  in  (  >  p^  (sG^^s)  )i) 


s  coefficient  of  sn  in  f(sG^_^(s)) 

Therefore,  Gk(s)  »  .  qk  nsn  -  f(sGk=1(s))0  Clearly,  gn(s)  a  sGn(s) 


Hi  0 


since  P(yn  =  n)  a  P(Sn  ~  n^  1)  and  gn(s)  -  sf(gna=1(s))0  Taking  the 
limit  as  n  approaches  infinity,  we  see  that  g(s)  3  sf(g(s))J 

Using  theorem  20,  the  expected  total  number  of  members  of  a 


family  is:  :?  (»)-  JVsfl! 
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and  the  variance  of  S  isi 
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ESTIMATION 


In  the  introduction  we  mentioned  that  the  parameters,  pr,  must 
be  empiraca  .ly  determined  for  the  given  amount  and  arrangement  of 
atomic  material.  The  method  of  maximum  likelihood  can  be  applied; 
however,  it  seems  necessary  to  assume  that  certain  random  variables 
are  observable. 

Let  be  the  number  of  individuals  in  the  m™1  generation 
who  have  exactly  r  descendents  in  the  (m4l)st  generation.  Then 

Let  S, 


oo 


— 


•Mr 


r*o 
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1  *1 


r=o  r-o 

Theorem  21.  The  maximum  likelihood  estimates  of  pr  and  JIX  based  on 
observed  values  of  z ^  for  min  (n^l  observed  generations)  are 

Ss  SL  7 _  Sn-f-i  ”  I 


Pr  = 


Zrrr 

Sn 


s, 


Proof:  In  order  to  construct  the  likelihood  function,  first  we 


determine  the  conditional  distribution  of  r  r  0,1, 


given 


“m 


From  our  initial  assumptions,  each  member  of  the  generation 


has  the  same  probability  distribution  associated  with  his  progeny* 
Therefore  the  probability  of  kQ  members  having  no  descendents  is 


/  Zm\  I  Zm 

P.  1  k  1  "here 


is  the  binomial  coefficient  and  is  equal  to 


•rr\ 


— yr--  .  Of  the  remaining  zm  -  kQ  members,  say  zm  -  k-^  have 
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one  descendent.  This  probability  must  be 
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Continuing 


in  this  manner  we  get  the  multinomial  functions: 


^ ( zmo>  zml> 


>  zmi> 


I  zm-l  ,0>  Vl.l 

oO 

7  I  TL  p 


« •  z 
9  m- 


l,i»  •**) 


77  {Zmr)\ 

Y-0 


“31- 


Therefore  P(znQ,  z^,  .  .  .  |  zn-l,0>  zn-l,l>  *  *  •>p(Bn-l,0>  Vl,l 


•  ' |zn-2,0  >  zn-2,l  » 
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=  TT  »  ,  .. 

"1  =  0  7 T  ( Z„  , ) ; 
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Due  to  the  Markov  properties  of  the  branching  process  this  is  the 
joint  distribution  of  the  z^.  ,  m  =  0,1,.  .  .  ,  r=  0,1,  .  ,  . 

The  only  factor  of  the  likelihood  function  which  depends  on  the  pr 
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TT  p..^- 7mr 


W3# 


oo  n  v 

the  logarithm  of  this  expression  is  ^  ; 

fro  M-0 


oo 


To  maximize  this  expression  subject  to  the  restraint,  ^  s | 


K=° 


the  method  of  Lagrangian  multipliers  may  be  used. 
™  L-Z  lzr»rl°S  P,  -Apr) 
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Setting  -  0  ,  we  get  vj  2wt,  - J pr  for  r  =  °>1s 
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But 


^  p  ~  0  5  summing  the  above  exoression  over  r  gives: 
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The  double  sum  of  the  is  merely  Sn  so  )  =  Sn  and 
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Using  these  facts,  it  follows  that: 
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ION  FLUCTUATION 


The  number  of  ion  pairs  produced  by  a  fast  primary  particle 
passing  through  an  absorbing  medium  is  a  branching  stochastic  process 
(Moyal  £l8]  ).  In  this  process  the  primary  particle,  which  may  be 
a  gamma  ray,  has  a  certain  charge,  mass,  and  velocity  or  a  given 
frequency;  the  thickness  and  atomic  properties  of  the  absorbing 
medium  are  known.  The  primary  interacts  with  the  atoms  in  the 
absorber,  freeing  electrons  by  the  Compton  effect0  These  electrons 
are  the  first  generation  and  each  may  ionize  another  atom  or  atoms 
depending  on  their  energy;  secondary  electrons  may  ionize  more 
atoms,  etc.  The  process  ends  when  all  the  freed  electrons  are  too 
slow  to  produce  further  ionization.  This  results  in  some  number  of 
free  electrons  and,  thereby,  an  equal  number  of  ion  pairs. 

In  determining  the  distribution  of  this  number  we  shall  make 
the  following  assumptions: 

1.  Successive  ionizing  collisions  are  statistically  independent. 

2.  The  total  energy  loss  of  the  primary  is  very  much  smaller 
than  its  initial  energy  and,  hence,  the  loss  of  primary  energy  may 
be  neglected.  Therefore,  we  are  not  treating  slow  primaries  or 
very  thick  absorbers.  We  also  assume  that  energy  losses  due  to 
radiation,  nuclear  interactions,  etc.  are  negligible. 

3.  The  absorbing  medium  is  homogeneous. 

In  The  recording  of  data  is  delayed  sufficiently  to  allow 
termination  of  the  process.  Photoelectric  effects  in  the  gas  or 
chamber  walls  is  neglected. 

Unlike  the  branching  processes  treated  in  the  previous  sections 


the  probability  distribution  of  the  number  of  progeny  per  individual 

in  a  given  generation  is  dependent  upon  the  generation,,  Thus  one 

may  expect  different  methods  to  be  used0  Moyal  has  applied  the 

Laplace  transform  method  to  solve  the  Kolmogorov  equations  of  a 

Markov  process  and  this  is  outlined  below. 

Let  t  be  the  thickness  of  the  absorber 

N  =  the  number  of  absorber  atoms  per  unit  volume 

ff(E)  -  total  cross  section 

q  =  primary  ionizing  collision  rate 

qn  =  P(n  ion  pairs  are  produced  in  a  given  collision) 

pft)  -  P(  n  ion  pairs  produced  in  a  thickness  t) 
oo 

M(s,t)  =  pn(t)e“ns  ,  the  generating  function  of  pn(t) 

n-o 

Then  (n  d’(E)  (dp)  dt)  is  the  probability  of  an  energy  loss  between 
E  qnd  G  -|-  dE  in  an  absorber  of  thickness  dt„  The  primary  ionization 
rate  is  c\  =  N  d£ 

OQ 

The  C|o  =  j  and  q(dt)  =  P(  an  ionizing  collision  in  thickness  dt) . 
n-o 

Since  at  least  one  ion  pair  is  produced  in  a  given  collision,  q0  =  0« 
Based  upon  the  above  assumptions  we  have 

<A>  Pnlt+ti)  =i; 

pn  ( gt)  =(I-Cj (St)  ^no  +  +  where  is  the 

Kronecker  delta  and  0(St)  goes  to  zero  faster  than  £t  as  £t*^0. 

By  assumption  1  and  equation  (A), 
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Comparing  this  equation  with  4f  -  =  £Y  >  11  is  clear  that 

M(*,t)=  exP[cit^(e‘"s-i)%  7  -  e®R(s) 

/)*(  J 

where  q  =  c,t 


n~i 


The  inversion  of  the  Laplace  transform,  M(s,t),  is  given  by: 


Pn  (Q)  ~  2 ttI  S 


C  +  TT  i* 


QRM 


no 
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c-m 


(see  Widder  [21]) 


Evaluating  this  integral  by  the  saddlepoint  method  (see  Korn  [l5]), 
using  the  first  term  only,  gives: 

P„(Q)  =  -g-[air  0  R’'(s„)]'/a e*p[d( Rls„)-  S„  R  (S„)j 

where  Sn  is  related  to  n  by 

n=  -QR'(Sn)  =  Q^>()<e~l<inJlch<} 

l<*> 
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since 


R(sJ  =2  (e  -lk< 

K-  I 


R"(S.)  -  f  ,<^K 
K-i 


-S„K 


and  c  is  a  normalizing  constant. 

Since  Z  p„(«)  =  I 

(0=0 

OO  .  VsS,  I 

c  =  2  D^ttQ R"(S„)J  exp  [Q{R(s„)-SnR  (s„ 

n-o 
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